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Abstract. We discuss an extension of the Hamilton- Jacobi theory to nonholonomic mechanics 
with a particular interest in its application to exactly integrating the equations of motion. We give 
an intrinsic proof of a nonholonomic analogue of the Hamilton-Jacobi theorem. Our intrinsic proof 
clarifies the difference from the conventional Hamilton-Jacobi theory for unconstrained systems. 
The proof also helps us identify a geometric meaning of the conditions on the solutions of the 
Hamilton-Jacobi equation that arise from nonholonomic constraints. The major advantage of our 
result is that it provides us with a method of integrating the equations of motion just as the 
unconstrained Hamilton-Jacobi theory does. In particular, we build on the work by Iglesias-Ponte, 
de Leon, and Martin de Diego [15] so that the conventional method of separation of variables applies 
to some nonholonomic mechanical systems. We also show a way to apply our result to systems to 
which separation of variables does not apply. 

1. Introduction 

1.1. The Hamilton-Jacobi Theory. The Hamilton-Jacobi theory for unconstrained systems is 
well understood both from the classical and geometric points of view. Besides its fundamental 
aspects, such as its relation to the action integral and generating functions of symplectic maps, the 
theory is known to be very useful in exactly integrating Hamilton's equations using the technique 
of separation of variables [see, e.g., [21 [HI [20]. See also Abraham and Marsden [1] Chapter 5] for 
an elegant geometric treatment of the Hamilton-Jacobi theory. 

1.2. Extension to Nonholonomic Mechanics. Our objective is to extend the Hamilton-Jacobi 
theory to nonholonomic systems, that is, mechanical systems with non-integrable velocity con- 
straints, building on the previous work by Iglesias-Ponte et al. p^. Nonholonomic mechanics deals 
with such systems by extending the ideas of Lagrangian and Hamiltonian mechanics [see, e.g., [5]. 
However it is often not straightforward to extend the ideas of unconstrained dynamics to nonholo- 
nomic systems, since a mechanical system loses some properties that are common to (conventional) 
Lagrangian and Hamiltonian systems when one adds nonholonomic constraints. 

Since the Hamilton-Jacobi theory is developed based on the Hamiltonian picture of dynamics, 
a natural starting point in extending the Hamilton-Jacobi theory to nonholonomic systems is a 
Hamiltonian formulation of nonholonomic mechanics. Bates and Sniatycki [3 and van der Schaft 
and Maschke [26 generalized the definition of Hamiltonian system to the almost-symplectic and 
almost-Poisson formulations, respectively [see also[5l [T71 [18]. As is shown in these papers, adding 
nonholonomic constraints to a Hamiltonian system renders the flow of the system non-symplectic. 
In fact, van der Schaft and Maschke [26 showed that the condition for the almost-Poisson Hamil- 
tonian system to be (strictly) Poisson is equivalent to the system being holonomic. This implies 
that the conventional Hamilton-Jacobi theory does not directly apply to nonholonomic mechanics, 
since the (strict) symplecticity is critical in the theory. In fact, the Hamilton-Jacobi equation is a 
PDE for generating functions that yield symplectic maps for the flows of the dynamics. 

There are some previous attempts to extend the Hamilton-Jacobi theory to nonholonomic me- 
chanics, such as Pavon [23]. However, as pointed out by Iglesias-Ponte et al. [15], these results are 
based on a variational approach, which does not apply to nonholonomic setting. See de Leon et al. 
[12j for details. 
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Iglesias-Ponte et al. [15] proved a nonholonomic Hamilton-Jacobi theorem that shares the geo- 
metric view with the unconstrained theory by Abraham and Marsden [T]. The recent work by 
de Leon et al. [12 developed a new geometric framework for systems defined with linear almost 
Poisson structures. Their result generalizes the Hamilton-Jacobi theory to the linear almost Poisson 
settings, and also specializes and provides geometric insights into nonholonomic mechanics. 

1.3. Nonholonomic Hamilton— Jacobi Theory. The previous work by Iglesias-Ponte et al. |15j 
and de Leon et al. [12] is of theoretical importance in its own right. However, it is still unknown 
if the theorems are applicable to the problem of exactly integrating the equations of motion of 
nonholonomic systems in a similar way to the conventional theory. To see this let us briefly discuss 
the difference between the unconstrained Hamilton-Jacobi equation and the nonholonomic ones 
mentioned above. First recall the conventional unconstrained theory: Let Q be a configuration 
space, T*Q be its cotangent bundle, and if : r*Q — > M be the Hamiltonian; then the Hamilton- 
Jacobi equation can be written as a single equation: 

//(,,^)=E, (la) 

or 

HodW{q)=E, (lb) 

for an unknown function W : Q ^ M.. On the other hand, the nonholonomic Hamilton-Jacobi 
equations in [15^ have the following form: 

d(Fo7)(g)GP°, (2) 

where 7 : Q ^ T*(5 is an unknown one- form, and T)° is the annihilator of the distribution T) C TQ 
defined by the nonholonomic constraints. While it is clear that Eq. ([2]) reduces to Eq. ([T]) for 
the special case that there are no constraint^ Eq. (j2]) in general gives a set of partial differential 
equations for 7 as opposed to a single equation like Eq. ([T]) . 

Having this difference in mind, let us now consider the following question: Is separation of 
variables applicable to the nonholonomic Hamilton-Jacobi equation? First recall how separation 
of variables works in the conventional setting: One first assumes that the function W can be split 
into pieces, each of which depends only on some subset of the variables q, e.g., 

W{q) = Wi{qi) + W2{q2), 

for Wi, W2 : Q — > M, and q = (gi, (72)- Then this sometimes helps us split the left-hand side of the 
Hamilton-Jacobi equation ([T]): 

with some functions Hi, H2 : T*Q — > M, and hence both Hi and H2 must be constant: 

where Ei and E2 are constants such that E1+E2 = E. Then we can solve them to obtain dWi/dqi 
and dW2/dq2 separately. It is not clear how this approach applies to the nonholonomic Hamilton- 
Jacobi Equation Furthermore, there are additional conditions on the solution 7 which do not 
exist in the conventional theory. 



V = TQ and hence 0° = and identifying the one-form 7 with dW 
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1.4. Integrability of Nonholonomic Systems. Integrability of Hamiltonian systems is an in- 
teresting question that has a close link with the Hamilton-Jacobi theory. For unconstrained Hamil- 
tonian systems, the Arnold-Liouville theorem [see, e.g., [2] stands as the definitive work. The link 
between the theorem and the Hamilton-Jacobi theory lies in the action-angle variables, which spec- 
ify the natural canonical coordinates for the invariant tori of the system; in practice the action-angle 
variables can be found through separation of variables for the Hamilton-Jacobi equation [see, e.g., 

m §6.2]. 

For nonholonomic mechanics, however, the Arnold-Liouville theorem does not directly apply, 
since the nonholonomic flow is not Hamiltonian and so the key ideas in the Arnold-Liouville theorem 
lose their effectiveness. Kozlov [19j gave certain conditions for integrability of nonholonomic systems 
with invariant measures. However, it is important to remark that there are examples that do not 
have invariant measures but are still integrable, such as the Chaplygin sleigh [see, e.g.,IH|5]. Also 
it is unknown how this result may be related to the nonholonomic Hamilton-Jacobi theory, which 
does not have an apparent relationship with invariant measures. 



1.5. Main Results. The goal of the present paper is to fill the gap between the unconstrained 
and nonholonomic Hamilton-Jacobi theory by showing applicability of separation of variables to 
nonholonomic systems, and also to discuss integrability of them. For that purpose, we would 
like to first reformulate the nonholonomic Hamilton-Jacobi theorem from an intrinsic point of 
vie"w|^ We show that the nonholonomic Hamilton-Jacobi equation ^ reduces to a single equation 
Hoj = E. This result resolves the differences between unconstrained and nonholonomic Hamilton- 
Jacobi equations mentioned in Section 1.3, and makes it possible to apply separation of variables 
to nonholonomic systems. Furthermore, the intrinsic proof helps us identify the difference from 
the unconstrained theory by Abraham and Marsden [1] and find the conditions on the solution 
7 arising from nonholonomic constraints that are more practical than (although equivalent to, as 
pointed out by Sosa [24 ) those of Iglesias-Ponte et al. [15]. It turns out that these conditions are 
not only useful in finding the solutions of the Hamilton-Jacobi equation by separation of variables, 
but also provide a way to integrate the equations of motion of a system to which separation of 
variables does not apply. 



1.6. Outline of the Paper. In Section [2] we briefiy review the Hamiltonian formulation of non- 
holonomic mechanics, and also state some definitions and results that pertain to the nonholonomic 
Hamilton-Jacobi theorem. In particular, we first give an intrinsic description of nonholonomic 
Hamilton equations, define and state a few results concerning completely nonholonomic constraints 
and regularity of nonholonomic systems. Much of the ideas in the proof of the nonholonomic 
Hamilton-Jacobi theorem come from identifying both the similarities and differences between the 
nonholonomic and unconstrained Hamilton equations. 

In Section |3] we formulate and prove the nonholonomic Hamilton-Jacobi theorem. The theorem 
and proof are an extension of the one by Abraham and Marsden [1 to the nonholonomic setting. In 
doing so we identify the differences from the unconstrained theory; this in turn gives the additional 
conditions arising from the nonholonomic constraints. 

We apply the nonholonomic Hamilton-Jacobi theorem to several examples in Section |4j We first 
apply the technique of separation of variables to solve the nonholonomic Hamilton-Jacobi equation 
to obtain exact solutions of the motions of the vertical rolling disk and knife edge on an inclined 
plane. We then take the snakeboard and Chaplygin sleigh as examples to which separation of 
variables does not apply, and show another way of employing the nonholonomic Hamilton-Jacobi 
theorem to exactly integrate the equations of motion. The conclusion follows to suggest possible 
future work. 



'A coordinate-based proof is given in [15] 
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2. Hamiltonian Formulation of Nonholonomic Mechanics 

2.1. Hamilton's Equations for Nonholonomic Systems. Hamiltonian approaches to nonholo- 
nomic mechanical systems are developed by, for example, Bates and Sniatycki f3] and van der Schaft 
and Maschke ^26j . See also Koon and Marsden [171 HH] and Bloch |5j| . 

Consider a mechanical system on a differentiable manifold Q with Lagrangian L : TQ M. 
Suppose that the system has nonholonomic constraints given by the distribution 

V:={veTQ\ uj'{v) = Ay = 0, s = 1, . . . ,p} , (3) 

where are Lagrange multipliers and lo^ = Af dq^ are linearly independent non-exact one-forms on 
Q. Then the Lagrange-d'Alembert principle gives the equation of motion [see, e.g. |5| Chapter 5]: 

d dL dL 
dtW-W = ^'^- 

The Legendre transformation of this set of equations gives the Hamiltonian formulation of non- 
holonomic systems. Specifically, define the Legendre transform FL : TQ T*Q hy 

d^ 
de 



¥L{Vq) ■ Wq = -j-L{Vq + £ Wq] 



e=0 

for Vq,Wq G TqQ. Throughout the paper we assume that the Lagrangian is hyperregular, i.e., the 
Legendre transform FL is a diffeomorphism. Set p := FL(g), or locally pi = dL/dq^, and define the 
Hamiltonian H :T*Q ^Rhy 

H{q,p) := {p,q) - L{q,q), 
where q = {¥L)^^{p) on the right-hand side. Then we can rewrite Eq. Q as follows: 

with the constraint equations 

'^'(9) = (^) =0 s = l,...,p. (6) 

Equations ^ and ^ define Hamilton's equations for nonholonomic systems. We can also write 
this system in the intrinsic form in the following way: Suppose that = q^dqi + Pi^p^ is the 
vector field on T*Q that defines the flow of the system, $7 is the standard symplectic form on T*Q, 
and ttq : T*Q ^ Q is the cotangent bundle projection. Then we can write Hamilton's equations 
for nonholonomic systems ([s]) and ^ in the following intrinsic form: 

ij^nhi^ = dH — XsTTqLo^ , (7) 

along with 

T7TQ{X^')eV or LO%T7TQ{Xf}'))=0ioi s = l,...,p. (8) 
Introducing the constrained momentum space M. := ¥L{'D) C T*Q, the above constraints may be 
replaced by the following: 

p£M. (9) 

2.2. Completely Nonholonomic Constraints. Let us introduce a special class of nonholonomic 
constraints that is assumed in the nonholonomic Hamilton-Jacobi theorerrEl 

Definition 2.1 (Vershik and Gershkovich [27/, see also Montgomery [[22]). A distribution D C TQ 
is said to be completely nonholonomic (or bracket- generating) if T> along with all of its iterated Lie 
brackets [T>,T>], [T>, [D,!)]], . . . spans the tangent bundle TQ. 

Let us also introduce the following notion for convenience: 
We would like to thank the referees for pointing out the importance of the notion. 
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Definition 2.2. Let Q be the configuration manifold of a mechanical system. Then nonholonomic 
constraints on the system are said to be completely nonholonomic if the distribution V C TQ 
defined by the nonholonomic constraints is completely nonholonomic (or bracket-generating). 

One of the most important results concerning completely nonholonomic distributions is the 
following 1^ 

Theorem 2.3 (Chow's Theorem). Let Q he a connected differentiahle manifold. If a distribution 
T) C TQ is completely nonholonomic, then any two points on Q can be joined by a horizontal path. 

We will need the following result that easily follows from Chow's Theorem: 

Proposition 2.4. Let Q be a connected differentiable manifold and D C TQ be a completely 
nonholonomic distribution. Then there is no non-zero exact one-form in the annihilator T>° C T*Q. 

Proof. Chow's Theorem says that, for any two points go and qi in Q, there exists a curve c : [0, T] 
Q with some T > such that c(0) = qo and c(T) = qi, and also c{t) £ ^c(t) for i ^ (0> T). Now 
let df be an exact one-form in the annihilator T>°. Then by Stokes' theorem, we have 

f{qi)-f{qo)= r df{c{t))dt = Q, 







where df{c{t)) = because df £ V° and c{t) G ^c(t)- Since go and qi are arbitrary and Q is 
connected, this implies that / is constant on Q. □ 

2.3. Regularity of Nonholonomic Systems. We will also need to assume regularity of non- 
holonomic systems in the following sens^ Consider a nonholonomic system with a hyperregular 
Lagrangian L : TQ M and a constant-dimensional distribution T> C TQ defined by nonholonomic 
constraints. For any Vq G TQ define a bilinear form BL{vq) : TqQ x TqQ ^ M by 

BL{Vg){Uq,Wq) := 7: ^ L{Vq -\- ElUq -\- £2Wq) 

0£iO£2 



= D2D2L{q,v) ■ {Uq,Wq). 

£l=£2=0 



Then hyperregularity of the Lagrangian implies that the associated map B^^{vq) : TqQ — > T*Q 
defined by 

{B^L{Vq){Uq),Wq'^ := B L{Vq) {Uq , Wq) 

is an isomorphism. Thus we can define a bilinear form Wl : T*Q x T*Q ^ M. hy 

WL{Vq)iaq,(3q) .= (oq , {Bi)-\(3q)] 



Definition 2.5 (de Leon and Martin de Diego [TT]; see also de Leon et al. [13]). In the above setup, 
suppose that the annihilator T>° is spanned by the one-forms {lo'^}'^^^. Then the nonholonomic 
system is said to be regular if the matrices (C£'*(u)) defined by 

Crfiv) = -WL{v){u;^,u;') (10) 

are nonsingular for any v £ T>. 

For a mechanical system whose Lagrangian is kinetic minus potential energy, regularity follows 
automatically: 

Proposition 2.6 (Carihena and Rahada |8]; see also de Leon and Martin de Diego [llj. // the 
Lagrangian L : TQ M has the form 

L{q,v) = ^g{v,v)-V{q), (11) 
with g being a Riemannian metric on Q, then the nonholonomic system is regular. 



^See, e.g., Montgomery 22 for a proof. 

^We again would like to thank one of the referees for pointing out the necessity of this assumption. 
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Proof. In this case I?2-D2L(g, v){ug, Wg) = gq{uq, Wq), and so Wl is defined by the inverse g"^^ of the 
matrix gij. Since gij is positive-definite, so is the inverse g^^; hence it follows that Wl is positive- 
definite. A positive-definite matrix restricted to a subspace is again positive-definite, and so is 
positive-definite and hence nondegenerate. □ 



In the Hamiltonian setting with the form of Lagrangian in Eq. ( 11 ), we have the following result: 



Proposition 2.7 (Bates and Sniatycki [3J. Suppose that the Lagrangian is of the form in Eq. (11). 
Let T he the distribution on T*Q defined by 

J^:={ve TT*Q I TTTQiv) G V} , (12) 

and then define a distribution Ti on M. := ¥L{V) by 

n:=J^nTM. (13) 

Then the standard symplectic form i7 restricted to TC is nondegenerate. 

Proof. See Bates and Sniatycki |i3j Theorem on p. 105]. □ 

3. NONHOLONOMIC HAMILTON-JACOBI THEOREM 



We would like to refine the result of Iglesias-Ponte et al. [15] with a particular attention to 
applications to exact integration of the equations of motion. Specifically, we would like to take 
an intrinsic approach (see [15j for the coordinate-based approach) to clarify the difference from 
the (unconstrained) Hamilton-Jacobi theorem of Abraham and Marsden [1] (Theorem 5.2.4). A 
significant difference from the result by Iglesias-Ponte et al. jl5j is that the nonholonomic Hamilton- 
Jacobi equation is given as a single algebraic equation H o ^ = E just as in the unconstrained 
Hamilton-Jacobi theory, as opposed to a set of differential equations d{H 07) . 

Theorem 3.1 (Nonholonomic Hamilton-Jacobi). Consider a nonholonomic system defined on 



a connected differentiable manifold Q with a Lagrangian of the form Eq. (11) and a completely 
nonholonomic constraint distribution D C TQ. Let ^ : Q ^ T*Q be a one- form that satisfies 

7(g) e Mq for any qeQ, (14) 

and 

c^tI'Dx'D = 0, i.e., d'y{v,w) = for any v,w £ v. (15) 
Then the following are equivalent: 

(i) For every curve c(t) in Q satisfying 

c{t) = T^Q-XH{loc{t)), (16) 

the curve 1 1— > ^oc{t) is an integral curve of X'^, where Xu is the Hamiltonian vector field 
of the unconstrained system with the same Hamiltonian, i.e., ix^^ — dH. 

(ii) The one-form 7 satisfies the nonholonomic Hamilton-Jacobi equation.- 

Hoj = E, (17) 

where E is a constant. 

The following lemma, which is a slight modification of Lemma 5.2.5 of Abraham and Marsden 
[1], is the key to the proof of the above theorem: 



Lemma 3.2. For any one-form j on Q that satisfies the condition Eq. (15) and any v,w £ J-, the 
following equality holds: 

i7(r(7 o vtq) -v.w) = 9.{v, w - r(7 o ttq) • w). (18) 
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Proof. Notice first that v — T(7 o ttq) • v is vertical for any v £ TT*Q: 

TlTQ ■ {V - r(7 O TTq) ■ v) = T-Kq{v) - T{'KQ 070 vtq) • V 

= Tttq{v) - Tttq{v) = 0, 
wliere we used the relation ttq 070 vtq = ttq. Hence 

9.{v - r(7 o ttq) • f , w - r(7 o ttq) • u;) = 0, 

and thus 

r2(r(7 o ttq) ■v,w) = Q{v, w - r(7 o hq) ■ w) + n(r(7 o ttq) ■ V, r(7 o ttq) • w). 

However, the second term on the right-hand side vanishes: 

r2(T(7 o ttq) • v,T{-f o ttq) • w) = ^*9.{T-kq{v),Tttq{w)) = -d'y{T-KQ{v) ,T-kq{w)) = 0, 

where we used the fact that for any one- form (5 on Q, = —dl3 with (3 on the left-hand side 
being regarded as a map (3 : Q ^ T*Q [SeelH Proposition 3.2.11 on p. 179], and the assumption 
that d'ylxyxV = 0; note that VjW £ implies T'Kq{v),T-kq{w) G P. □ 

Let us state another lemma: 

Lemma 3.3. The unconstrained Hamiltonian vector field Xh evaluated on the constrained mo- 
mentum space M. is in the distribution T , i.e., 

Xnioig) £ Tocq for any Og G Mg. 

Proof. We want to show that TTrQ(XH{aq)) is in Vg. First notice that 

OH 

TTTQ{XH{aq)) = Q^iaq) = ^H{aq), 
where we defined ¥H ■.T*Q ^TQ by 

d 
de 

However, because the Lagrangian L is hyperregular, we have ¥H = (¥L)~^ and thus 

TTTQiXHiaq)) = (¥Lr\aq). 

Now, by the definition of Ai, aq £ M implies ag £ ¥L(T>g), which gives (FL)-i(ag) £ Vg by the 
hyperregularity of L. Hence the claim follows. □ 

Proof of Theorem 3.1 Let us first show that implies Assume and let p{t) := 7 o c(t), 
where c{t) satisfies Eq. (16). Then 

pit) = T^{c{t)) 

= T-f ■ Tttq ■ Xni'y o c{t)) 

= T{jOTTQ)-XH{7 0c{t)). (19) 

Therefore, using Lemmas |3.2| and |3.3[ we obtain, for any w £ 

0(r(7 o ttq) • XH{p{t)),w) = n{XH{p{t)),w - r(7 o ng) ■ w) 

= n{XH{p{t)),w)-n{XH{p{t)),T{^0 7TQ)-w). 

For the first term on the right-hand side, notice that for any w £ J^, 

niX'^'^w) = dH-w-XsTT*Qio'{w) = dH ■w = Q{Xh,w). 
Also for the second term, 

n{XH{p{t)),T{^ o ttq) • w) = dH{p{t)) ■ r(7 o ttq) • u; = d{H o ^){c{t)) ■ Tttq{w). 



{Pg,¥H{ag)) = H{ag+£Pg] 



e=0 
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So we now have 

n{T{j o ttq) ■ XH{p{t)),w) = n{Xf{p{t)),w) - d{H o 7)(c(t)) • TTTQiw). (20) 

However, the nonholonomic Hamilton- J acobi equation ( |17[ ) imphes that the second term on the 
right-hand side vanishes. Thus we have 

n{T{j0 7TQ)-XH{p{t)),w) 



(21) 



for any w £ ^p{t)- Now T(707rQ) ■ Xh G TM. since 7 takes values in A^; also T{'yoTTQ) ■ Xnipit)) G 



J^p[t) because 

TvTQ • r(7 o ttq) • XH{p{t)) = TiTTQ o 7 o vtq) • Xnipit)) = Tttq ■ Xnipit)) G V, 

using Lemma 3.3 again. Therefore T(7 o ttq) ■ Xnipit)) G Hpi^ty On the other hand, X^'^(p(t)) G 
'Hp(^t) s-s well: Xjj^{p{t)) G Tp^^-^M because M. is an invariant manifold of the nonholonomic flow 
defined by X'f and also X^^(p(t)) G Tp^t) due to Eq. (g. Now, in Eq. ([2T|, w is an arbitrary 



element in Tp(t) and thus Eq. \2l\ holds for any u; G 'Hp(t) because H d T . However, according to 
Proposition 2.7 restricted to H is nondegenerate. So we obtain 



r(7 o ttq) • Xu{p{t)) 



x^nh 



and hence Eq. ( 19 ) gives 



p{t)=Xf{p{t)). 

This means that p{t) gives an integral curve of X^. Thus ^ implies (p]). 

Conversely, assume let c{t) be a curve in Q that satisfies Eq. (16) and set p{t) := 7 o 
Then p{t) is an integral curve of Xfj^ and so 



cit). 



Pit) = xftm) 



However, from the definition of p(t) and Eq. ( |16| ), 

p(t) = r7(c(t)) = r7 • Tttq • Xnipit)) = o vtq) • X^^(p(t)). 

Therefore we get 

xf{p{t)) = n^oT,Q)-XHm)- 



In view of Eq. (20), we get, for any w G TT*Q such that Tirqiw) G P, 

(i(//o7)(c(t))-r^Q(u;) =0, 

but this implies d{H o 7)(c(t)) • u = for any f G I?c(t)) oi' '^(-f^ ° 7)(c(t)) G ^°(()- However, this 
further implies d{H o j){q) = for any q £ Q: For an arbitrary point q £ Q, consider a curve c{t) 
that satisfies Eq. ( |16[ ) such that c(0) = Then this gives d{H o^)[q) £ 'D°. Therefore d{Ho^) £ T>° 
on Q, but then Proposition 2.4 implies that d{H o 7) = because T) is assumed to be completely 
nonholonomic. Therefore we have H o ^ = E for some constant E, which is the nonholonomic 
Hamilton- Jacobi equation (17). □ 



Remark 3.4. The condition on ^7, Eq. (15), stated in the above theorem is equivalent to the one 
in [15j as pointed out by Sosa [2lj [see also |22l Lemma 4.6 on p. 51]. However Eq. (15) gives a 



simpler geometric interpretation and also is easily implemented in applications. To be specific, the 
condition in [T5] states that there exist one-forms such that 

d7 = X]rAa;^ (22) 

s=l 

which does not easily translate into direct expressions for the conditions on 7. On the other hand, 

d'y{vi,Vj) = for any i / j, (23) 



Eq. (15) is equivalent to 
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where {vi}^Zi spans the distribution T>. Clearly the above equations give direct expressions for the 
conditions on 7. We will see later in Section [4] that the above equations play an important role in 
exact integration. 

Remark 3.5. Table [T] compares Theorem |3.1| with the unconstrained Hamilton- Jacobi theorem 
of Abraham and Marsden [1^ (Theorem 5.2.4). Note that Eq. (15) is trivially satisfied for the 
unconstrained case. Recall that 7 is replaced by an exact one-form dW in this case. Since D = TQ 
by assumption, we have d^\x>xv = dj = d{dW) = and thus this does not impose any condition 
on dW. 

Table 1. Comparison between unconstrained and nonholonomic Hamilton- Jacobi theorems. 





Unconstrained 


Nonholonomic 


Generating Function 


W -.Q^R 


None 


One-form 


dW -.Q^ T*Q 


-.Q^ M CT*Q 


Condition 


ddW = (trivial) 


d-ylvxv = 


Hamilton-Jacobi Eq. 


HodW{q) = E 


Hoj{q) = E 



Remark 3.6. See Carihena et al. [S] for a Lagrangian version of Theorem 3.1 and de Leon et al. 
P^2j for an extension to a more general framework, i.e., systems defined with linear almost Poisson 
structures. 



4. Application to Exactly Integrating Equations of Motion 
4.1. Applying the Nonholonomic Hamilton Jacobi Theorem to Exact Integration. The- 



orem 



3.1 suggests a way to use the solution of the Hamilton-Jacobi equation to integrate the 



equations of motion. Namely, 

Step 1. Find a solution 7(5) of the Hamilton-Jacobi equation 

Hoj{q) = E, 

that satisfies the conditions 7(g) G M.q and dj\x>xv = 0; 



(24) 



Step 2. Substitute the solution 7(g) into Eq. (16) to obtain the set of first-order ODEs defined in 
the configuration Q: 

c{t) = TTTQ-XHi-focit)), (25a) 

or, in coordinates. 



c{t) 



dH 

dp 



(7oc(t)); 



(25b) 



Step 3. Solve the ODEs (25) to find the curve c(t) in the configuration space Q. Then 7oc(t) gives 
the dynamics in the phase space T*Q. 

Figure [l] depicts the idea of this procedure. 

In the following sections, we apply this procedure to several examples of nonholonomic systems. 
In any of the examples to follow, it is easy to check that the constraints are completely nonholo- 



nomic (see Definition 2.1), and also that the Lagrangian takes the form in Eq. (11) and hence the 



system is regular in the sense of Definition 2.5 
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Figure 1 . Schematic of an implication of the nonholonomic Hamilton-Jacobi theorem. 

z 

A 




Figure 2. Vertical rolling disk. 



4.2. Examples with Separation of Variables. Let us first illustrate through a very simple 
example how the above procedure works with the method of separation of variables. 

Example 4.1 (The vertical rolling disk). [See, e.g., |5]. Consider the motion of the vertical rolling 
disk of radius R shown in Fig. [2} The configuration space is Q = SE{2) x = {(x, y, 99, •i/')}- 
Suppose that m is the mass of the disk, / is the moment of inertia of the disk about the axis 
perpendicular to the plane of the disk, and J is the moment of inertia about an axis in the plane 
of the disk (both axes passing through the disk's center). The velocity constraints are 

X = R cos ifTp, y = R sin f Tp, (26) 

or in terms of constraint one-forms, 

uj"^ = dx — R cos (fdil), uj'^ = dy — RsiiKf dip. (27) 

The Hamiltonian H : T*Q is given by 

2 \ m J II 
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The nonholonomic Hamilton- Jacobi equation (17) is 



(29) 



where is a constant (the total energy). Let us construct an ansatz for Eq. (29). The momentum 
constraint p ^ ^A gives px = mi? cos ipp^/ 1 and py = mR sin ipp^/I, and so we can write 7 : Q ^ 
M as 



mR 



mR 



7 = — cos If j^{x,y, If, '4))dx + — sin (p-f^{x,y,ip, 4)) dy 



+ Ifix, y, if, V') dip + 7^(x, y, if, ijj) dijj (30) 



Now we assume the following ansatz: 



Then the condition d'yl'pxv = in Eq. (15) gives 

dip 



0, 



and so 



So Eq. (29) becomes 



l^p{x,y,(p,'ip) = 7^(x,y,V')- 
-7^(a;,y, V-) 



E. 



(31) 

(32) 
(33) 
(34) 



2 V J ' P 

The first term in the parentheses depends only on ip, whereas the second depends on x, y, and tp. 
This implies that both of them must be constant: 

= 7°, y, V') = 7° , (35) 

where 7^ and 7^ are the constants determined by the initial condition such that 

1 / 1 , n.o / + mR'^ 



E. 



Then Eq. ( 16 ) becomes 



X = , COS 



"Or 



y = sm 99, 



/ " I 

which are integrated easily to give the solution 



x{t) = Cl + 

y{t) = C2 



7° 
J 



(36) 



JRil 


sin 


(1% 


It 




JRil 


cos 


(1 






7° 
J ' 


m = 



(37) 



4 
I 



where ci, C2, V'o, and ^/^o are all constants. 



Separation of variables for unconstrained Hamilton-Jacobi equations often deals with problems 
with potential forces, e.g., a harmonic oscillator and the Kepler problem. Let us show that separa- 
tion of variables works also for the following simple nonholonomic system with a potential force: 
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Figure 3. Knife edge on inclined plane. 

Example 4.2 (The Knife Edge; see, e.g., Bloch [5]). Consider a plane slanted at an angle a from 
the horizontal and let (x, y) represent the position of the point of contact of the knife edge with 
respect to a fixed Cartesian coordinate system on the plane (see Fig.[3|. The configuration space is 
Q = SE(2) = {{x, y, (p)}. Suppose that the mass of the knife edge is m, and the moment of inertia 
about the axis perpendicular to the inclined plane through its contact point is J. The velocity 
constraint is 

simpx — cosipy = 0, (38) 

and so the constraint one-form is 

uj^ = sin Lf dx — cos (f dy . (39) 

The Hamiltonian H : T*Q ^ R is given by 

\ I -\- \ 
^ _ _ I _x_ y_ _|_ _^ 1 _ ^gxsiua. (40) 



The nonholonomic Hamilton- Jacobi equation (17) is 



Hoj = E, (41) 



where ii^ is a constant (the total energy). Let us construct an ansatz for Eq. (41 ). The momentum 
constraint p £ M gives 

Py = tan ifp^, 

and so we can write j : Q ^ M as 

7 = Jx{x,y,ip)dx + tanip-/cc{x,y,ip)dy + 'j^{x,y,ip)dip. (42) 
Now we assume the following ansatz: 

j^{x,y,<f) =-f^{(f). (43) 
Then the condition d'ylxixv = in Eq. (15) gives 

^ = -tan(/9 7^. (44) 

Integration of this equation yields 

Ixix, y, ^) = f{x, y) cos V9, (45) 
with some function f{x,y). Then Eq. (|41[) becomes 



1 



{2mgsina) x + 



m J 



E. (46) 
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The first two terms in tlie brackets depend only on x and y, whereas the third depends only on ip. 
This implies that 

= 7° (47) 

with some constant 7^, and f{x,y) satisfies 



fix,yr ,^ . ^ All? 
- [2mg sm a) x + 



m 



J 



E. 



(48) 



Let us suppose that sleigh is sliding downward in Fig. [3] Then we should have ^ for < 99 < 
7r/2. From Eq. (45) we see that f{x,y) > 0, and hence choose the branch 



f{x,y) 



m 2E 



J 



+ (2m'^gsma) x. 



(49) 



Then Eq. ( 16 ) becomes 



cos if 




(7^ 
2J 



+ {mg sin a) x, 



sm If 
V = — , 

ym72 

Let us choose the initial condition 




(50) 



(7°)- 
2J 



+ (mg sin a) 



hp 
J ' 



(x(0),y(0),v9(0),i(0),y(0),vi(0)) = (0, 0, 0, 0, 0, u;), 



where uj := 7^/ J- Then we obtain 
xit) 



sm [ujt), y{t) 



2tj2 



gfsina / 1 , 
2u;2 ( ~ ^sin(2a;t) ) , ip{t) = ut. 



(51) 



These are the solution obtained in Bloch [5l Section 1.6]. 



4.3. Examples without Separation of Variables. In the unconstrained theory, separation of 
variables seems to be the only practical way of solving the Hamilton-Jacobi equation. However 
notice that separation of variables implies the existence of conserved quantities (or at least one) 
independent of the Hamiltonian, which often turn out to be the momentum maps arising from 
the symmetry of the system. This means that the integrability argument based on separation of 
variables is possible only if there are sufficient number of conserved quantities independent of the 
Hamiltonian [See, e.g.,|20l §VIII.3]. This is consistent with the Arnold-Liouville theorem, and as a 
matter of fact, separation of variables can be used to identify the action-angle variables [See, e.g., 
Ea §6.2]. 

The above two examples show that we have a similar situation on the nonholonomic side as well. 
In each of these two examples we found conserved quantities (which are not the Hamiltonian) from 
the Hamilton-Jacobi equation by separation of variables as in the unconstrained theory. So again 
the existence of sufficient number of conserved quantities is necessary for application of separation 
variables. However, this condition can be more restrictive for nonholonomic systems since, for 
nonholonomic systems, momentum maps are replaced by momentum equations, which in general 
do not give conservation laws [6 . 

An interesting question to ask is then: What can we do when separation of variables does not 
seem to be working? In the unconstrained theory, there are cases where one can come up with a 
new set of coordinates in which one can apply separation of variables. An example is the use of 
elliptic coordinates in the problem of attraction by two fixed centers [21 §47. C]. The question of 
existence of such coordinates for nonholonomic examples is interesting to consider. However, we 
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would like to take a different approach based on what we already have. Namely we illustrate how 
the nonholonomic Hamilton-Jacobi theorem can be used for those examples to which we cannot 
apply separation of variables. The key idea is to utilize the condition d'^\vy.v = 0, which does not 
exist in the unconstrained theory as shown in Remark |3.5[ 

Example 4.3 (The Snakeboard; see, e.g., Bloch et al. Pi). Consider the motion of the snakeboard 
shown in Fig. |4j Let m be the total mass of the board, J the inertia of the board, Jo the inertia of 




Figure 4. The Snakeboard. 

the rotor, Ji the inertia of each of the wheels, and assume the relation J + Jq + 2 Ji = mr^. The 
configuration space is Q = SE(2) x x = {{x, y, 9, jp, ^)} and the Hamiltonian H : T*Q ^ U. 
is given by 

The velocity constraints are 

X + r cot (j) COS0 6 = 0, y + r cot (p sinO = 0, (53) 
and thus the constraint distribution is written as 

p = = (x, y, 9, ij, (/>) G TQ I uj%v) = 0, s = 1, 2} , (54) 

where 

uj^ = dx + r cot cj) cos 9 d9, liP' = dy + r cot (j) sin 9 d9. (55) 
The nonholonomic Hamilton- Jacobi equation (17) is 

Hoj = E. (56) 



Let us construct an ansatz for Eq. (56). The momentum constraint p ^ Ai gives 

TYIT TTIT 

Px = n T cot (p COS 9 {pe - p^), py = 5 cot (j) sin 9 {jpe - p^)-, 

mr-^ — Jo V- » — Jo 

and so we can write j : Q ^ Ai as 

TflT 

7 = ^ — cot (j) hg - 7^,) (cos 9dx + sin 9 dy) + jgd9 + 7^, dip + 7^ d(p (57) 

mr^ — Jo 

Now we assume the following ansatz: 

li>{x,y,9,Tp,cl)) ='y^{ip), -f^{x,y,9,'tp,(p) = -/^{(j)). (58) 



Then the nonholonomic Hamilton-Jacobi equation (56) becomes 

2(mr2 - Jo)2 '^(^^-TV) + + 2(mr2 - Jo) " ^u[^^ = ^- 

Solving this for 70, we have 



(59) 



^e{x,y,e,iP,^)=-1^{'^) + —=======^E-^ ^— (60) 

V(mr^ - Jo sin^(/))/2 V ^"^o 4Ji 
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and substituting the result and Eq. (58) into the condition (i7|DxX' = in Eq. (15) gives 

d 



sm ( 



Therefore it follows that 



I'M = 0- 



for some constants 7^ and 7^. Hence Eq. (4.3) becomes 

(mr^ — Jo) C sin (j) 



where we defined 



C 



'E 



0^2 



Then Eq. ( 16 ) gives 





2 Jo 4Ji 

C r cos cos ( 



5(0) 

g((/)) := \J {mr"^ - Josin^ 0)/2 



C r sin ^ cos < 



5('A) 



C sin ( 



7^ C sin ( 



2Ji' 



(61) 



9(0) ' ffC*/*) 

This result is consistent with that of Koon and Marsden obtained by reduction of Hamilton's 
equations for nonholonomic systems. It is also clear from the above expressions that the solution 
is obtained by a quadrature. 

In the above example we found conserved quantities through d'ylxyxv = instead of separation of 
variables. In the following example, we cannot identify conserved quantities even through djlxixv = 
0; nevertheless we can still integrate the equations of motion. 

Example 4.4 (The Chaplygin sleigh; see, e.g., Bloch [5 ). Consider the motion of the Chaplygin 
sleigh shown in Fig. [5} Let m be the mass, / the moment of inertia about the center of mass C, a 




Figure 5. The Chaplygin sleigh. 
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be the distance from the center of mass C to the contact point A of the edge. The configuration 
space \s Q = SE{2) = {(x, y, 0)}, where the coordinates (x, y) give the position of the contact point 
of the edge (not the center of mass). The velocity constraint is 



sin 9 X — cos 9y = 0, 



and so the constraint one-form is 



The Hamiltonian H :T*Q 



LJ^ = sin 9 dx — cos 9 dy. 
is given by 



(62) 
(63) 



Ma^ sin^ 9 + J o Ma^ cos^ 9 + J o 
H = pi + WTT7 Pv 



2JM 



2JM 



1 

2J' 



a sm 9 cos 9 <^ , n 

J PxPy + -j{sm. 9 - cos 9 py)pe. (64) 



The nonholonomic Hamilton- Jacobi equation (17) is 



(65) 



where ii^ is a constant (the total energy). Let us construct an ansatz for Eq. (65). The momentum 
constraint p ^ M. gives 

aM sec 9 



Py = tant'px + 



J + a?M 



P9, 



and so we can write j : Q ^ Ai as 

7 = 7x{x,y,9)dx + 
Now we assume the following ansatz: 



aM sec 9 

tan 9 -fx{x, y, 9) + — — ^ 7e(x, y, 9) 



leix,y,9) = -fe{9). 



dy + jg{x,y,9)d9. 



Then the condition d^\vxv = in Eq. (15) gives 



{J + a^M)sec9 ( + tan 6*73; ) + aM tan 6^ ( +tan6l70 



dje 



d9 



0. 



(66) 
(67) 
(68) 



On the other hand, the Hamilton-Jacobi equation (65) becomes 
"2 sec (9 



1 



sec 9 



M 



, ^,9 4a tan 6' 
7^(x,y,6' +— — 27-.lx{x,y,9)-ie{9) 
J + a^M 



(J + 2a^M + J cos 29) sec 9 
(J + a2M)2 



7e(^)' 



(69) 



It is impossible to separate the variables as we did in the examples in Examples 4.1 and 4.3 



since 



we cannot isolate the terms that depend only on 9. Instead we solve the above equation for 7^; and 
substitute the result into Eq. (68). Then we obtain 



-aWM 2E 



Solving this ODE gives 



drie_ 
d9 



7g(6') = (J + a^M) cjcos 



J + a?M ' 



am 

J + a^M 



(70) 
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where we assumed that x'{0) = y'{0) = 0, 6{0) = 0, and ^'(0) = lo and also that \9{t)\ < tt/2; note 
that the angular velocity lo is related to the total energy by the equation E = [J + a'^M) uj/2. Then 



the equation for 9{t) in Eq. (16) becomes 



= i^coslx -^^^e] , (71) 



J + a2M 



which, with 6{0) = 0, gives 



2 

6(t) = - arctan 
b 



tanh - ujt 



\2 

where we set b := \J c^Mji^J + a'^M). Substituting this back into Eq. (70), we obtain 



(72) 



Je{t) = {J + a^M)uJsedl(^^JJ^^u;tj, (73) 

which is the solution obtained by Bloch [3] [see also|5l Section 8.6]. 

5. Conclusion and Future Work 

We formulated a nonholonomic Hamilton-Jacobi theorem building on the work by Iglesias-Ponte 
et al. |T5] with a particular interest in the application to exactly integrating the equations of motion 
of nonholonomic mechanical systems. In particular we formulated the theorem so that the technique 
of separation of variables applies as in the unconstrained theory. We illustrated how this works 
for the vertical rolling disk and snakeboard. Furthermore, we proposed another way of exactly 
integrating the equations of motion without using separation of variables. 

The following topics are interesting to consider for future work: 

• Relation between measure-preservation and applicability of separations of variables. The in- 
tegr ability conditions of nonholonomic systems formulated by Kozlov [19j include measure- 
preservation. As mentioned above, applicability of separation of variables implies the ex- 
istence of conserved quantities other than the Hamiltonian. Therefore it is interesting to 
see how these ideas, i.e., measure- preservation, applicability of separation of variables, and 
existence of conserved quantities, are related to each other. 

• "Right" coordinates in nonholonomic Hamilton- J acobi theory and relation to quasivelocities. 
In the unconstrained Hamilton-Jacobi theory, there are examples which are solvable by 
separation of variables only after a certain coordinate transformation. As a matter of fact, 
Lanczos \2Gi\ p. 243] says "The separable nature of a problem constitutes no inherent feature 
of the physical properties of a mechanical system, but is entirely a matter of the right system 
of coordinates^ It is reasonable to expect the same situation in nonholonomic Hamilton- 
Jacobi theory. In fact the equations of nonholonomic mechanics take simpler forms with 
the quasivelocities [3 [10]. Relating the "right" coordinates, if any, to the quasivelocities is 
an interesting question to consider. 

• Extension to Dirac mechanics. Implicit Lagrangian/Hamiltonian systems defined with 
Dirac structures |25| l28l |29] can incorporate more general constraints than nonholonomic 
constraints including those from degenerate Lagrangians and Hamiltonians, and give non- 
holonomic mechanics as a special case. A generalization of the Hamilton-Jacobi theory to 
such systems is in progress [21]. 
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